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Formulas are developed for the diffracted field around a perfectly conducting cylinder 
with a dielectric sleeve of arbitrary thickness. These formulas represent the field due to a 
unit electric line source parallel to the cylinder (either inside or outside the dielectric sleeve) 
as a spectrum of radial eigenf unctions. It is shown that in each case the field in the rc^gion 
containing the source can t)e expressed as the sum of two terms, the first of which is a spectral 
representation of the field w^hen the outer as well as the inner surface of the dielectric is a 
perfectly conducting boimdary. The second term of the sum, which alone involves the 
properties of the dielectric, is an integral that converges rapidly at high frequencies. Using 
these general results, perturbation calculations are carried out for three limiting examples 
of plane wave scattering: (1) Thin dielectric; (2) low-refractive index; (3) very small sur- 
face curvature. In the latt(M- n])])roxiination the correspondence with o])tical results is 
shown. 

1. Introduction 

The scattei-ing of j)lano oloclTomagnetic wavos by cylincHcal obstacles has received detailed 
theoretical treaiment only in the two special cases for wliicli the scalar wave function, ^, 
satisfies either tlie Dirichlet or the Neumann boundary condition at the cylindrical surface. 
In electromagnetic terms these problems corres])()n(l rc^spectively to the scattering by a per- 
fectly conducting cyliinU'r of a wave whose electi-ic or juagnetic vector is parallel to the cylinder 
axis. For a cylindrical obstacle wliose siu'face is not perfectly conducting a jnore complicated 
boundary condition of tlu^. type ^ + ad,,^ = is recjuired (sometimes called an impedance 
boundary conchtion). in such easels the theory has not l)een so extensively elaborated.^ 

Til is paper is concerned with some of the thc^oretical features of the scattering problem 
associated with a cylindrical obslach' of this more gcMu^'al kin(L in ])articular we shall examine 
tlie scattered wave procUiced by a (Ueleclric cylinckT of circular cross section having a perfectly 
conducting coaxial core. One aspect of this problcMu has Ijcen investigated experimentally by 
Tang [1] ^ for an incident wave with electric vector polarized parallel to the cylinder, and in 
the same paper the Foiuicn' expansion of the w^ave ftmction has been used for numerical calcu- 
lations. It seems, however, that no general discussion of tiie characteristics of the wave 
solution has been given; it is this task which will now occupy oin* attention. 

2. Derivation of the Green's Function 

The cylinder which we shall consider (fig. 1) is composed of a perfectly conchicting- core 
of radius a surrounded by a dielectric layer of inner radius a and outer radius h. In the sub- 
sequent analysis the relative dielectric constant of the layer will be taken to be a real scalar 
gi'eater than one, er=€/eo>l. That is to say, the refractive index of the dielectric is a real 
number greater than that which characterizes the surrounding medium (air or vacuum) . For 
a lossy or anisotropic dielectric the theory is unaffected except by the substitution throughout 
of a complex or tensor dielectric constant for e. Of course, the detailed calculation of field 
quantities will be much more complicated in such cases. 

If the composite cylinder described above is excited by a harmonic line current parallel 
to its axis (electric line source, exp[— twt]), the task of finding the resultant electric field at 
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Figure 1. The composite cyl- 
inder and coordinate system. 



every point outside the perfect conductor resolves itself into the determination of the Green's 
functions which satisfy appropriate partial differential equations and boundary conditions. 
For the moment we shall limit our considerations to the case in which the line source is outside 
the dielectric. Denoting coplanar position vectors relative to an origin on the cylinder axis 
by r' and r for source and observation points, respectively, we require solutions in the xy-iplsme 
for the scalar wave equations 



(V'+F)ffi(r,r0=0; a<r<b,r'>b 
subject to the following boundary conditions at r=a and r=b: 

^i(a,rO-0; 
^o(b,rO = ^i(b,rO; 



(2.1a) 
(2.1b) 

(2.2a) 
(2.2b) 
(2.2c) 



Here the wave number, k=co/v, is related to the dielectric constant and the refractive index, 
n, of the dielectric according to 

k=ko\er=nko. (2.3) 

The mathematical problem defined by eqs (2.1) and (2.2) may be solved in a variety of 
ways. We shall employ a technique developed and discussed by Titchmarsh [2] and Marcu- 
vitz [3] for separable equations. The procedure consists in constructing so-called characteristic 
functions which satisfy the related system of ordinary differential equations 






e-]' 



\e',X) = -8{d-d'). 



(2.4a) 
(2.4b) 
(2.4c) 



These equations are found from (2.1) by separating the polar variables {r, 6) with X taking the 
role of separation constant. In addition to (2.4) the functions /o and/i are required to satisfy 
boundary conditions corresponding to (2.2), while the angular function g is made unique by 
certain convergence conditions to be specified. Once the characteristic functions have been 
explicitly determined from (2.4) and the boundary conditions, the Green's functions required 
by (2.1) are obtained by an integration in the complex X-plane. For it can be verified that if 
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C is a circular contour of iniiiiito radius in this plane, then 

^•^^dXf(ry,\)=r8(r-n, (2.5a) 

^. j)jl\g {d,d^\)=b (e-e^) . (2.5b) 

With the help of (2.4) and (2.5) direct substitution in (2.1) will sliow that the required Green's 
function is 

Gir,t')=^.j^d\f(r,r',\)cj(e,e',X), (2.6) 

where F is a contour which encloses the singularities of one of the functions /or g, the sense of 
the contour being given either by (2.5a) or (2.5b) as the case may be. 

In carrying out the procedure outlined above, we shall begin by considering the charac- 
teristic function which satisfies (2.4c). It has the form 

ry(^,^^X)=.46-(^-V, (2.7) 

where v'^=X and dy, 0< are respectively tlie larger and the smallei' of the two variables 6 and d\ 
The positive exponential solution has been chosen to ensure the regularity of g for all values of 
X on the principal sheet of the two-sheeted X surface. On this slieet Tm ^^^O, and since 

0< arg p= arg VX <7r, (2.8) 

the region of regularity corresponds to the upper half of the ^-plane. The constant A in (2.7) 
is determined by the fact that dg/dd has unit discontinuity at 6=6^ for a delta function soui'ce. 
We thus find from (2.4c) that the characteristic function is 

^(^,^^X)=^6-^-V = ^^^''''-''^ < arg . <7r. (2.9) 

It is easy to verify that (2.9) satisfies (2.5b). Since the oidy singularity of g(d, d\ X) is a 
branch point at X=0, the clockwise contour Cmay be deformed to an integral around the branch 
cut which is conveniently taken along the positive real axis of the X-plane. If the deformed 
contour is now mapped onto the real axis of the j^-plane by the transformation X=v^, the 
resulting integral is seen to be a representation of tlie 5-f unction. Thus 



^.^./x^(^,^^x)=^J_\/.6^>'^-^''=5(^-o. 



The characteristic Green's functions which satisfy (2.4a) and (2.4b) are constructed in a 
similar fashion but from linear combinations of Bessel functions of order v=^J \ . At large dis- 
tances from the source, r:>r\ they must represent outgoing waves. Furthermore, the interior 
function must vanish at the conducting surface r = a, and the exterior function must behave 
like a ^'standing wave" in the region r^y>r>b. In terms of Hankel functions of the first and 
second kinds the most general expressions wliicli exhibit these characteristics are 

Ury,\) = B,[m'^ {hh)H?' {hr<)-DH?^ {hh)m' {hr<d]Hl'^ (hr>); ry>h, (2.10) 

Mr/ ,\) = B,[m'' (ka)m'' (kr^) -H?' {ka)TIi'' {kr^Wl'' (k^ry) ; a<r< h, r' >b. (2.11) 

Two of the three integration constants appearing in (2.10) and (2.11) are determined by the 
pair of linear algebraic equations which express the conthiuity of the radial function and its 

first derivative at r= 6. The remaining constant is fixed by the jump in r-r-fo(;ry,\) at r=r\ 
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Making use of the Wi-onskian relation 

i?<» (^c)ff«' '(ft) -H«> ' (ft)iy<» (ft) =-^ , (2.12) 

and introducing the notation 

a=ka, ^=kb, Po=hb; (2.13) 

'C,(kr) =m'^ {a)Hi'' (kr) -7?J2) (^)^(i) (/.^) ^ 

C^^ytr) ^fl^i^^ (a)Zf ^2) / (^^) _£f (2) (^);;^(i) ^ (^^) . (2.14) 

and 

V(1.2)_/o ^'-W , /3 ^L_^^Li§K, (9 1^^ 

we find that the characteristic function for an exterior source at r' is 

fo{r,r',X)='^^H'/\hr<)-^^^m'Khr<)jm'\kor>), r>b; (2.16) 

^'^"^'"^ '^^-wmcMW' - ■ ^^ ^ 

A direct calculation shows that these functions satisfy (2.5a). 

It remains to determine (?(r,r') from (2.6). First we note that the only singularities of 
/o and/i in the X-plane are simple poles at the set of complex values, \j, for which 

X^J?-0. (2.18) 

That these functions have no branch point is easy to verify. With the relations, 

H^}l=e'^''Hl'' and H^^l=e-'^''m'\ 

we learn that 

C-,= a and X_,=X. 

It follows that/o and /i are both single-valued over the entire X-plane. In accordance with 
(2.6), we now form the integral 

(?(r,r') =2^. J^^/(r,r',X) e''^<'-^\ (2.19) 

where the contour r is chosen to enclose the singularities of/ in a positive sense. Since the 
product function has a branch point at X=0, T must be indented around the branch cut, which 
we shall take along the real axis as shown in figure 2. In the j/-plane the resulting contour 
maps to 7. By virtue of (2.8) and the as^^mptotic behavior of the Hankel functions, the contri- 
bution to (2.19) from the semicircular segment of 7 may be made arbitrarily small if the radius 
of the arc is made lar^e enough. In this way we see that both a discrete and a continuous 
representation of G are possible, 

^(r.rO-S Kes (^=^j)=^f dvf{r/,\) e-'^-^'L (2.20) 



5 It would be more explicit to write C^ia, kr), but since subsequent operations on C^ affect it only as a function of kr, it is perhaps clearer to 
use the abbreviated notation. 
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u- plane 




Figure 2. The complex v (uui X planes. 



X- plane 



3. The Discrete Spectrum 

The discrete representation of the Green's function is obtained by evaluating the residues 
of f{ry,\) from (2.16) and (2.17). A more symmetrical form results if one makes use of tlie 
relation 

4i 1 



X!^''=Xl'' + 



which, at v=vj, becomes simply 



X?)_ 



4i 



TT 7/yUA.)//:^H/?o) 



1 



TiiriMHi'^m 



C^.i) 



(:^.2) 



wliere the subscript indicates the j^^ root of (2.18). With this change the residue sums corre- 
sponding to (2.19) form discrete spectral representations of the total field due to an exterior line 
source at r' > b. Thev are 






(^i (.r,r j -1 z^ ^^^<., CW)H^'' (/3o) ' ''-'-''• 



(3.3<i) 



(3.3b) 



It is now clear from coiisidsrations of symmetry that if the lino source is inside the dielectric, 
a<r' <b, than 






(3.4a) 
(3.4b) 



It is iiotcwortliy that the representations (3.3a) and (3.4b), with source and observation points 
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in the same region, are unaffected by an interchange of the variables r and r' . If r and r' are 
in different regions, however, their interchange turns (3.3b) into (3.4a) and vice versa. 

The radial eigenfunction expansions derived above have certain characteristic features 
which should be mentioned. For fixed r and r' each term of the series represents (in the 
terminology of Franz) a ^ ^creeping wave" propagating in the ^-direction with exponentially 
damped amphtude. Since d covers the interval (—0°, qq), the total field at an}- point is the 
sum of the contributions at ^di2m7r for all integral values of m. At high frequencies, however, 
the imaginary part of i^ is a large positive number so that all terms other than the one for which 
m=0 are exponentially small. Asymptotically, (3.3) and (3.4) represent the total field in the 
range 0<|^— 0'|<27r. It is to be expected, moreover, that the set of Vj are ordered with in- 
creasingly large imaginary parts so that if |^— ^'| ^0, only the first term of each residue expan- 
sion will contribute significantly to the field at high enough frequencies. However, the problem 
of finding even the first root of (2.18) would seem to be a formidable one. 

4. The Continuous Spectrum 

An alternative representation of the Greenes function as an integral over a continuous 
spectrum is obtained from (2.19), using (2.16) and (2.17). Since g{d^ d\ X) has a branch point 
singularity at X=0, the contour integral becomes 

G(T,r')=-^CJt^e^-^'-'y{ry,\), (4.1) 

where the indicated path is around a branch cut along the real axis of the X-plane. With the 
help of (3.1) and the change of variable, \=v^, we find from (4.1) that 

Go(r,r')=~ J^ dv[mH^^,)H?\hr<)-m'^\^,)m\hr<)] ^Jir^ e^'''''' 

Similarly the field inside the dielectric is 

^^^'''^-2xJ_A'i^> Cmm\M ' r>b,a<r<b. (4.3) 

In this form the Greenes function of (4.2) has some interesting features. The first term 
satisfies all the conditions associated with the field of an electric line source in the vicinity of 
a perfectly conducting cylinder of radius b when the medium is characterized by the propaga- 
tion number ko. One way to see this is to close the integration contour in the upper half of the 
j'-plane. The residue sum that results is 

^ rj(2) fo \ 

Gc(r,r')=^ S ^PtI) H'v\k^<)m\hr>)e'''''-'\ (4.4) 

where 

m'\M=mi'{M=o- (4-5) 

From the Wronskian, (2.12), one finds that 
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so Unit (4.4) can be rearranged in the form 

These are the results for a conducting cyHnder which have been derived in other ways by T. 
Wu [4] and by the author [5]. 

We see then that the field of a composite cylinder has as one part the creeping waves of 
a conducting cylinder. Only the second integral of (4.2) depends upon the parameters of the 
dielectric (thickness and refraction index) and this dependence resides in the factor Xj,^\ 
Indeed as a-^b, Xj,^^ grows without limit and the integral vanishes, leaving the Green's function 
for a conducting cylinder of radius b. On the other hand, as k-^ko the limiting value of the 
second integral of (4.2) is 

and on substituting we find that the Green's function reduces to that for a conducthig cylinder 
of radius a. These limiting forms will be used later for perturbation calculations. 

When tlie line source is hiside (lie dielectric sleeve, the derivation of tlie Green's function 
proceeds in a similar way. We begin with characteristic functions for the radial equations 
having tlie form 

h,{ry,\)=A,r,{kr<:)Hl'\h;ry)', a</ <b, r>b, (4.9) 

h^{r/,\)^Ar(\{kr^)\TI?\kry) + imi'\kry)]] a<ry<b. (4.10) 

Both functions satisfy tlu^ boiiudary coiidilioii at r<:^ = a, and hp represents an outgoing cylin- 
drical wave as r>— >oo . The constants can be evaluated as before, and after a certain amount 
of manipulation, we find that 

/, (rr' M 1 (\{kr')m\hr) . . 

h,{r/,\) = ~'^i^-^[W^\m?{kr>)-n?\mi'\h'>)]+^ ^^ (4.12) 

The exterior function hpiry ,\) of (4.11) corresponds exactly to/i of (2.17) after an interchange 
of the variables r and r', so the related Green's function corresponds to (4.3). When both 
source and observation points are inside the dielectric, the convolution of hi(ry ,\) with (2.9) 
around the branch cut of g(dy,\) leads to the following integral representation of the Green's 
function : 

\ n dv C.{kr<)CXkr>) ,,,,_,., ,.,.,. 

The first term of this Green's function vanishes at r<=a and ry = b, and one can easily verify 
that it corresponds exactl}^ to the field of an electric line source inside a dielectric filled coaxial 
transmission line with perfectly conducting walls. The second term represents the modification 
of that field due to radiation into the exterior medium. 
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The residue sum associated with (4.11) has already been given in (3.4). To complete the 
discussion we shall write down the symmetrical representations for a coaxial line, analogous to 
(4.7) and (4.8). They follow from the convolution of the first term of (4.12) with g(d,6\\) 
around the simple poles located at the points of the X-plane for which 

It is found that 

where we have made use of the Wronskian relations, 

m^'^ (kr) a (kr) -m^'' ' {kr) ft (^^) = ~^^^'''^ W • (4-16) 

5. Plane- Wave Excitation 

In order to render more manageable the integrals appearing in (4.2) it is convenient at 
this stage to turn one's attention to special cases. Since the general class of problems asso- 
ciated with the scattering of a plane wave is of considerable interest, we shall first consider the 
question of reducing the field to that of an incident plane wave of unit amplitude. This may 
be done by considering the first term of (4.2), which we shall call Gjy, 

(?.(r,r') =^ J"^ dv [i7<« (^or<) -gpjil H«> (Ar„r<)] H<'> (/c„r>)<?'>i^-<"i. (5.1) 

The Hankel functions, Hj,{p), grow as exp( i^ In — j when j/-^+°° , so the individual terms of this 

integral diverge and cannot be treated separate^. However the divergences may be removed 
by first putting (5.1) in the form, 

G,(r,r^) =~ J^^ dv \h?^ (kor<) -§w^ ^'^ (^or<)] ^'^ (hr^) cos .(^-^0 , (5.2) 

and then making the substitution, 

^«>(p) = 2j,(p)-7?<»(p). 
With these changes it is found that 

G,{r,T')=^-j^^dvMhr<:)m'Hkor>) cos ^i^-^')-\f^d^4^) Hl'\hr<)m'\hr>) cos v{e-e'). 

(5.3) 
Each of these integrals is absolutely convergent. 

To establish a correspondence between (5.3) and the required plane wave solution, we 
follow Wu's procedure and sum the creeping waves for the first term. Then by interchanging 
the order of summation and integration (or by invoking the Poisson summation formula) it 
can be shown that 

\ iZ r dvJ,(kor<)m'^(kory) cos v{d-d' + 2m7r)==lm'\h\ry-T^\), (5.4) 

^ w=-co Jo 4 
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wlicrc tlie function on the right is just the free space Green's function representing a unit hne 
s()ui(*(\ If we now let r>— >oo and make use of tlie asymptotic form, 



H'^\p)' 



[2 i 

-■\—e 



(-l) 



both sides of (5.4) can be normahzed to an incident phme wave of unit ampUtude by multi- 
plying with the factor— 4'i(2)"^^^(/:o^>)exp—'i!(/:o^>—x/4). Making these changes in (5.3) and 
setting ^^ = gives the total field, i/'6(r), when a plane wave from the right is scattered by a con- 
ducting cylinder of radius b. This field is 

Mr) = 2 r dvMkor)e~''^ cos ve-2 f"' dv 4j^ m'^hr) e~ ''~^ cos vd. (5,5) 

Jo Jo ^p \Po) 

The procedure that led to (5.4) will also show that 

m= — co Jo 

Wlien (4.2) has been normalized in accordance with (5.5), tlie total field outside a dielectric 
clad cylinder illuminated by a plane wave from (lu^ riglit becomes 

,M /M^2 r- dv m'\hr) u{e-^) 
^(r)-Wr)+^J_^a)|HFTAO?' ' ^''^^ 

The last term is the modification of (lie scattered field introduced by the dielectric. It can also 
be put hito the form, 

which converges for all values of 6. 

6. The Surface Distribution 

In this section the distribution of ip on the dielectric surface will be examined in a number of 
special cases. Tliis choice of field point does not affect the methods of approximation that will 
be used, but does simplify the general expression for x// to a certain extent. For on the circle 
r=b the first term of (5.7) vanishes, so that the integral to be considered is simply 



^(•^)=^r 



x^'m^'\M 



e<'-^\ (6.1) 



In terms of the index of refraction, n=k/ko, and the quantity, 8 = 13— a, tlie conditions under 
which approximations to (6.1) will be investigated are (a) d/p<^l; (b) |?i— 1|<C1; and 
(c) koa:$>l, 

6.1. Thin Dielectric: 8/p<^l 

Since all the dielectric effects are contained in the function, 

we require the expansion of IjXy when a^b. It is clearly advantageous to treat Xy as a function 
of the variable a (i. e., the hmer radius a) since that parameter appears only in Cy. Moreover, 
H\}^ (Po) has no zeros for real v, so tlie first term of (6.2) predominates for small 8. We may 
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therefore put l/X^^Fp, where 



Since lim F^{a) = 0, the Taylor expansion of (6.3) has the form 

a->;S 

FXP-b) = -bF:{ei)+\ d'F'/i^)- .... (6.4) 

With (6.3) and (6.4) it is a straightforward matter to develop the required series exphcitly. 
We need certain relations of the Wronskian type [6], and up to terms of order 5^ the result of 
the calculation turns out to be^ 

i^.W=^+| 1+^1+^1+ • • • +ljAP'-^')+ljA2fi^-^r^)+ .... (6.5) 

It is clear that the first group of terms represents the function — ln[l — (5/i8)], so that 

On substituting (6.6) in (6.1) a series of integrals is obtained, the first of which is 

V a/TrO-ooi^^^H^o) 

The calculation of c)r\l/b(h) from (5.5) shows that except for the constant, In (6/a), the last result 
is just the current distribution on the surface of a perfectly conducting cylinder of radius b. 
Therefore, to a first approximation for thin dielectric sleeves of any refractive index and 
radius, the electric field distribution on the dielectric surface is proportional to the magnetic 
field distribution on a conducting cylinder of the same size when illuminated by a plane wave. 
The first term involving the refractive index is 0{8l^y and is proportional to n^. 

6.2. Low-Refractive Index: \n~l\<^l 
In order to expand the surface distribution, 

in the neighborhood of ko, let us consider the function, 

G.{k) = lXi'^m'\^o)]-'. (6.7) 

Its Taylor expansion may be written in the form 

GM = G.(ko) + in--l)koG:(ko) + 0(n-iy, (6.8) 

where {n-l)kQ=k—ko, Using (6.2) and the Wronskian (2.12), we find at once that 

^ . . __7^^ CJa,M _7ri I" „(2) fo \ Hj"^^ (o^o) ttcd ^r\'] {a q\ 

^^^^^^-T ^F(^-T L^^ ^^^^"FF^ ^. ^(^o)J- (6.9) 



' The successive derivatives of h\{a), evaluated at a=/3, are 
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Thus tlie first term in the expansion of (6.1) is just \paW, the total field at r—b when the plane 
wave is scattered by a conducting cylinder of radius a in the medium ko. 
The second term is found by evaluating 

where dA; signifies the differential operator d/dk. Treating Cy as a function of two dependent 
variables a, fi, and replacing CI by d^Cy, we have 

It follows, with the help of (6.9), that 

G'Xh) =~ ^§Sj^z [bC.^0C\+^C.d,di^,-0d,C\b/M^,^ (6.11) 

Since c)jc=ac>a^bc)^, the final bracket of (6.11) becomes, after a little rearrangement, 

I=kab{C.c)aC>^C-d^C,dar.)-ph{d^C;)'+bC.{bpC\+pdlC;), (6.12) 

If we now make use of BessePs differential equation for C„ the last parenthesis of (6.12) 
is seen to be 

d^a+/3d^a=^ {v'-^')a. (6.13) 

Moreover, the bracketed expression in the first term of (6.12) can be represented by the 
Wronskian determinant 

W{a,b/^)= . (6.14) 

In the usual way [7] one can show that 

which implies that pW is independent of 13. To find the value of the constant we let P = a, 
Then C,=0, and 

iim/3ir=-a(d«ad^a)/3=a. 

In the limit the derivatives are ordinary Wronskians like (2.12), and it follows that 

W=-^- (6.15) 

The result of substituting (6.13) and (6.15) into (6.12) is 

With this result, eq (6.8), (6.9), and (6.11) can be combined into an explicit expression for the 
first term of the expansion of G,(k) near k = ko. Substitution in (6.1) yields the following 
perturbed distribution at r=6 when n^l: 

+i£«-''>[^^]"'""*)«'-<'"^'<''}- <«■•« 

29 



If one is interested in the perturbation introduced by the dielectric at points outside the 
sleeve, it is only necessary to put Hi,^\kor) for H^^\Po) wherever that function appears explicitly 
in (6.16) . Wliether the field point is on the surface or outside the sleeve, however, it is readily 
verifiable that all the integrals converge and that the bracket vanishes when b = a, leaving 
only the field of the conducting cylinder. In general, it seems that none of the correction 
terms in (6.16) can be neglected. 

6.3. Very-High Frequencies: koCt^ 1 

The last method of approximation to be dealt with is useful when botli inner and outer 
cylinders are large compared to the wavelength in either medium. By combining the longer 
wavelength with the shorter radius these conditions are satisfied by the single relation, ao^l. 
When this requirement is met, no further restrictions need be imposed, a priori, on either the 
thickness of the dielectric sleeve or on its refractive index. Since we shall consider only the 
contribution to (6.1) from that part of the integral for which |j^! + 3|^|^^^<q:oj il is possible to 
use the Debye approximation for all Hankel functions. To this order of accuracy we may 
expect that the integration will yield the optical rays, but no diffraction effects. 

For a typical Hankel function the appropriate asymptotic form is 

m'^ (M -2'^'e-'kirPo sin yjo)-'^' exp ifioKvo), (6.17) 

where 770 is the real angle 

Vo=cos-'(v/fio), (6.18) 

and 

/()7o) -=sin r]o—Vo COS 77u. (6.19) 

The corresponding approximation for the derivative with respect to the argument, (So, is 

m'^\M-ismrjom'\M^ (6.20) 

With (6.17) and (6.20) it follows that 

P,-jf^j^iA sin vo. (6.21) 

In a similar way, by defining the angles 

r} = cos-'(v/p); ^=cos-'(p/a), (6.22) 

and the function 

F{v)=^f(v)~-af{t), (6.23) 

one can show that 

C'(3) .iF(p)_l -iF(p) 

If (6.21) and (6.24) are now combined according to (6.2), we find that 

1 1 ip l-e''""^'^ 



X /^osinT/ol+pl + re^^"'^''^ 
where p and r signify the ratios. 



(6.25) 



P=-^—' J T=T~\ — (6.26) 

13 sm rj l-\-p ' ^ 

Now for all j^^</5o, |S sin tj^ISq sin 770, so that 0>p>l . Consequently, 0<r<l and \t exp 2iF{v) |<1 . 
The denominator of (6.25) can therefore be represented by its Taylor expansion, and it follows 
that 

Tf«~«-^T^ [l-(l + ^)«"^'''+^l + ^)^^''^*' -^^ + ^)^'^'*''''''+- • •]■ (6.27) 
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Til is result already makes evident tlie ray luiture of the approximation. Successive 
terms of the series correspond to the rays which are, (1) reflected by the sleeve, (2) reflected 
once at the conducting cylinder, (3) reflected once at the conducting cyhnder and once inter- 
nally at the dielectric surface, etc. The phase function of each term corresponds to the number 
of path segments a ray traverses in the dielectric (with a minus sign added for each reflection 
at the conducting cylinder), and its amplitude function includes the appropriate transmission 
and reflection coefficients for each contact with the dielectric surface. Since these quantities 
are all functions of the integration variable, successive terms can be evaluated in detail only 
by more and more intricate stationary phase approximations. There is, however, one limiting 
case for which it is relatively easy to establish the correspondence between the terms of (6.27) 
and optical rays. When the radii a and h become indefinitely large, while their difference, 
d=b—a, remains constant, the solution for the cylindrical problem at each point h={b,d) of the 
dielectric surface approaches that for a plane wave scattered by a conducting plane sheet 
with a dielectric coating of thickness d, the dielectric surface being tangent to the cylinder at 
the point in question. 

To see how this comes about let us substitute (6.27) and (6.17) in (6.1). The first integral 
of the resulting asymptotic series becomes 

(6.28) 
the Debye approximations 

(6.29) 

For large /3o the principal contributions to this integral come from the neighborliood of the 
stationary points of g'-^^iv) in {—A, A). Since the asymptotic integrand decays exponentially 
outside this interval, one may extend its limits to (— «», «») without appreciable error. Then 
if a stationary point of g^^\p) is located at v=vi, tlie Taylor expansion of (6.28) in the vicinity 
of this point gives as its first approximation, 
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^^^ ^(^T ^"' ''"^'''" T+^ ''^J" ''" '^P ' [•^/"'(''')+^ ^^f/'^C".)]' (6.30) 



where \=v—vi and the subscripts indicate functions that are to be evaluated at v=vi. 

The value of v required in (6.30) is that whicli reduces c)^g^^\vi)=0 to an identity. It 
can be shown that this occurs for the angle (^70)1= 2^—^, corresponding to 

p, = fioSind (6.31) 

According to the condition specified prior to (6.17), the angle 6 must therefore be such that 

sin d+313-'/' (sin ey'<~ (6.32) 

It is clear that this inequality cannot under any circumstances be satisfied at the shadow 
boundary where d=7rl2. With (6.31) and (6.26) we determine that 

and 

l + p-^ocos^+0cos^, '^''^ ^^-^'^^ 

where R is the reflection coefficient for a plane dielectric surface [8] and the angles of incidence 
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and refraction are related by 

sin 6 



sin dj 
With these results, we have 



=n. 



^<" (b) ~i (1 +R) e-'^o cos . [^^J' Vf £ dX exp [- ,: ^^J (6.34) 

The value of the infinite integral is given in standard tables, and (6.34) can be reduced to 

iP^'^{h)^{l + R)e-'^o-ose^ (6.35) 



which corresponds in amplitude and phase to the sum of the incident and reflected rays at the 
dielectric interface for 
are no longer required. 



dielectric interface for angles ^<^ (fig. 3). Note that further restrictions on 0, such as (6.32), 




♦^ FiGUEE 3. Depicting ray geometry. 



The second term in the expansion of (6.1) is 

^«)(b) = _(^yVi J^^ci.(^„sin W-''^(j|^ exp ^fif«' W, (6-36) 

where 

^«'(.) = 2FW+s,"'W. (6.37) 

Again, after finding the point V2 at which (>yg^^\v2) vanishes, we can approximate (6.36) with 

^«'(b) (^y''(sin rjohjp^.e^zj^jX exp i^g''\v,)+^blg<'Kv,)j, (6.38) 

where \=v—v2 and the subscripts indicate functions that are to be evaluated at v=v2. 
According to (6.37) the required stationary point is the solution of 

c>A2F(v)+g''\v)]=0. (6.39) 

With the help of (6.23) and (6.29) it may be shown that this equation is equivalent to 

-2{n-^) + rjo^^T-d. (6.40) 

In order to discuss the solution of (6.40) in the plane limit we shall assume that ()? — f) is small so 
that sin (??— f) is approximately equal to its argument. Then 

r,-^^^^^. (6.41) 

sin ij /3 ^ ' 
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when' b = k(l and terms of order (5//3)^ have been neglected. Accordingly, if (ry — f) is to be a 
small quantity it is sufficient that d/p—d/b<^<^l, provided ryp^O or what amounts to the same 
thing, v'^9^0^. This requirement is, of course, the same as that already imposed by the original 
Debye approximations. With this understanding, (6.40) can be put into the more convenient 
form 

cos [770-2 ^?^^1=sin e, (6.42) 

L sm r? |8J 

Expanding the left-hand side of this relation up to terms of order bjfi, we obtain 



V2—vi= — 2piVi-j (6.43) 



wliere as before 

O ' n /3o COS i9 

Since in this limit V2 is a perturbation on vj, all functions of V2 in (6.38) have Taylor expansions 
in the neighborhood of vi. Retaining terms in 5//3 if tliey affect the phase of \p^^^ but not other- 
wise, we find, for example, that 

--^0 cos d+2{fi cos d,d ^+0 Qy. (6.44) 

To this degree of approximation (().3<S) })ecomes 

V'\h) \ (l + «)(l-/0 [ ^^,, COS J cxp[-i^„cos&+2i5cos&„l€'^^J_"^ d\ ...xp [-i^^^^), 

and we finally obtain (6.45) 

lA^'^b)'- (1+70(1-/?) oxp [-'/ft cos e+2ib cos 6>J, (6.46) 



where 



2^cos6>^ _ 

^-^-ftcos^+^cos^„-^^- ^^-^^^ 



The quantity, T^, is the transmission coefficient of a ray incident on the dielectric surface from 
the interior at an angle dn with the normal. If the transmission coefficient from the external 
region is designated by T=l + /?, then the sum of (6.35) and (6.46) has the form 

i/'(b)-6-^"^o^««n7^-77;e2^«^^^^"+ . . .]. (6.48) 

These two rays are depicted in fi,gure 3. It is an elementary exercise in geometry to show that 
their amplitudes and phases are exactly those predicted by ray optics, provided the dielectric 
surface is nearly plane over the arc segment illuminated by the bundle included between rays 
1 and 2. Additional rays for this quasi-plane example may be written down by inspection from 
the further terms of (6.27). 
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